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ABSTRACT 

This paper presents four searches for flaring sources of neutrinos using the IceCube 
neutrino telescope. For the first time, a search is performed over the entire parameter 
space of energy, direction and time with sensitivity to neutrino flares lasting between 
20 IIS and a year duration from astrophysical sources. Searches which integrate over 
time are less sensitive to flares because they are affected by a larger background of 
atmospheric neutrinos and muons that can be reduced by the use of additional timing 
information. Flaring sources considered here, such as active galactic nuclei, soft gamma- 
ray repeaters and gamma-ray bursts, are promising candidate neutrino emitters. 

Two searches are " untriggered" in the sense that they look for any possible flare in 
the entire sky and from a predefined catalog of sources from which photon fiares have 
been recorded. The other two searches are triggered by multi-wavelength information 
on flares from blazars and from a soft gamma-ray repeater. One triggered search uses 
lightcurves from Fermi-LAT which provides continuous monitoring. A second triggered 
search uses information where the flux states have been measured only for short periods 
of time near the flares. The untriggered searches use data taken by 40 strings of IceCube 
between Apr 5, 2008 and May 20, 2009. The triggered searches also use data taken by 
the 22-string configuration of IceCube operating between May 31, 2007 and Apr 5, 2008. 
The results from all four searches are compatible with a fluctuation of the background. 

Subject headings: triggered searches, multi-wavelength campaigns, blazars, soft-gamma 
ray repeaters 
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1. Introduction 

High-energy neutrinos can be produced in the direct vicinity of charged cosmic ray sources by 
the interaction of the high-energy cosmic rays with matter or photon fields. In those processes, 
charged pions are produced and result in a flux of high-energy neutrinos. The latter are unique 
messengers with which to observe the universe, as they have no charge and interact weakly, traveling 
directly from their point of creation essentially without absorption, distinguishing them from cosmic 
rays and high energy photons. Neutrinos are key in understanding the mechanisms of cosmic ray 
acceleration, and their detection from an astrophysical point source would be a clear indication 
of hadronic acceleration in that source. Time-integrated analyses suffer from a high irreducible 
background of atmospheric neutrinos and atmospheric muons, making them less sensitive for the 
detection of flares. Time-dependent analyses aim to reduce this background by searching over 
smaller time scales around the flares. The searches discussed in this paper are about a factor of 
four to five more powerful than time-integrated searches for flares of ~ 1 second. 

Four searches for time-dependent neutrino emissions from various categories of flaring sources 
are presented in this paper using data from the incomplete 22 and 40-string configurations of Ice- 
Cube. We call triggered searches those using multi-wavelength (MWL) information from photon 
experiments as a method of selecting flaring periods. Based on this information, we select cata- 
logues of interesting candidate neutrino-emitting flar es. We focus o n flares of larger durat i on tha n 



GRBs that are cove red by other IceCube searches lAbbasi et al.l (|2010bl ). lAbbasi et al.l (|2009c 



(jAbbasi et al.ll2011al )). Sources we consider are Active Galactic Nuclei (AGN) and one period of 



activity from the newly discovered Soft Gamma Repeater, SGR 0501+4516. 

The underlying assumption of triggered searches is that the neutrino emission follows the time- 
dependent emission of the photons, as a consequence of an excited state of the source when the jet 
can accelerate particles to higher energies than in the quiescent state. This hypothesis is assumed 
in the likelihood method as a prior. In order to make the search for flares as general as possible, an 
"untriggered" scan for clusters of events in time and direction without prior timing information is 
also performed. Untriggered searches are capable of detecting flares of duration similar to Gamma- 
Ray Bursts (GRBs) to flares lasting many days, as seen in AGNs. While dedicated GRB searches 
are more sensitive, the untriggered all-sky search is capable of detecting 'quenched' GRBs which 
may have been undetected by photon telescopes and hence not included in the dedicated searches. 
All-sky searches are affected by large trial factors, hence a catalogue of promising sources which 
are variable in MWL observations is selected for an additional untriggered search. 

The paper is structured as follows. In Section [2] the properties of the flaring sources, AGNs, 
Gamma- Ray Bursts (GRBs) and SGRs are considered. In Section [3] the data samples of 40 strings 
and 22 strings of IceCube are described. In Section |4] the time-dependent likelihood method is 
illustrated in general and compared with the time-integrated method. The four searches for flares 
are: 
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• An untriggered all-sky scan for short-duration neutrino emission from point-like sources (Sec- 
tion [5|); 

• An untriggered search for flares from a predefined catalogue of 40 sources identified as variable 
in GeV photons (Section [6]); 

• A triggered search for continuously monitored sources using MWL information. Data from 
the Fermi-LAT and also from SWIFT (Section [7]) are used; 

• A triggered search using sporadic information on flares collected by various X-ray experiments 
and Imaging Cherenkov Telescopes (lACTs) with typical energy thresholds around 300 GeV 
(Section [8|). These experiments have a reduced duty cycle compared to Fermi and tend to 
monitor sources mainly during their flaring states. 

In sections 5 to 8 each search method and its expected discovery potential are presented, and 
results are provided. In all searches the background is estimated directly from scrambled data, 
since the signal contribution is expected to be small. To avoid any bias toward discovery, each 
search has been performed in a blind fashion by defining cuts and search methods before looking at 
the final event sample. The final significance, which accounts for the different trials, is calculated 
by scrambling data in time. For each search the resulting probabilities (p-values) that the data 
could be due to background fluctuation are provided (see also Section U]). Trial factors between 
searches are not included in the p- values, as these trial factors are negligible compared to the 5a 
significance required by IceCube for a claim of discovery. 



2. Candidate Sources of Flaring Neutrino Emission 



Galactic and extra-galactic sources exhibit time-dependent variability ranging from short 
bursts with durations between seconds and minutes (e.g. GRBs or giant flares from SGRs) to 
long periods of high activity lasting hours to weeks (e.g. AGN flares). One of the main targets of 
the searches presented in this paper is variable emission from AGNs. Flat-Spectrum Radio Quasars 
(FSRQs) and BL Lacs, commonly unified in the AGN class of blazars, exhibit relativistic jets point- 
ing towards the Earth and some of the most violent variable high energy phenomena. Their spectral 
energy distributions (SEDs) extend orders of magnitudes across the electromagnetic spectrum and 
are characterized by structures of low and high energy non-thermal peaks. After 11 months o f 
operation, the Fermi-LAT collaboration published their first AGN catalogue (jAbdo et al.l l2010l ) 
containing 709 GeV-sources a ssociated with AGj Vs, many of which are in the previously published 
Bright Source list catalogue (lAbdo et al.l l2009bl^. A previous stacking search for neutrinos from 
AGNs used AMANDA data (jAchterberg et ahlboOGcl ). 



The low energy component in the radio to soft X-rays is due to synchrotron radiation of 
electrons gyrating in a magnetic field. The high energy component (X-ray to 7-ray) is explained 
in leptonic models by synchrotron emissions of electrons in the jet and subsequent up-scattering 
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of photons (inverse Compton) by the electron population responsible for the synchrotron emission. 
Hadronic models are used to explain the high energy component with relativistic protons, with 
energies above the threshold for p-7 or p-p pion production, which will decay to 7-rays and also 
neutrinos. Proton synchrotron emission can also contribute to the hi gh-energy component if the y 
are accelerated to ver y high energies (for a review on models see e.g. iNeronov &l Ribordvl (|2009l ): 
Bottched ((20101, 120071 ) and references therein) . 



The emission from blazars is known to be variable at all wavel engths. Simultaneous MWL ob- 
servations are crucial for understanding the cause of this variability (iGaidos et al. I ll99d :l iBlazejowski et al 
2OO5I : iKartaltepe k Balonekl I2OO7I : iHoran et aPbood : iBottcher et aD bood ) . The intensity of these 
objects can vary by more than an order of magnitude between different observing epochs. The 
typical time scales of AGN flares vary from hours to days, though high-energy variabili ty has been 
observed on much short er time scales, in some cases even down to just a few minutes (jAharonian 
2OO7I : I Albert et al.ll2007l ^. 



While leptonic models enjoy relatively good success reproducing the observed emission, other 
arguments favor a hadronic component. Perhaps the most compelling evidence are observations of 
"orphan" flares, defined as TeV photon eniission without accompanying X-rays, such as the lES 
1959-^650 flar e in 2002 (IKrawczvnski et ahlboO-d ). An a posteriori observation with AMANDA-II 



of two events (Bernardini et al 



200 



one exactly during the flare and ariother 31 days later, trig- 



2005 



Reimer et al 



2010 



gered some theoretical calculations (jHalzen &: Hooper 

included in the MW L triggered searches (see S ections [7] and [8] also (jCiprini 



2OO5I'). Two recent flares 



2005 



Abdo fc others 



Acciari et al.ll2009al : IVittorini et al.ll2009l )) are suspected to be orphan flares, but X-ray ob- 



servations were not simultaneous with gamma-ray observations and there is a possibility of having 
missed the X-ray flare. 



GRBs, believ ed to be produced by the most powerful phenomena in the uni verse (IMeszaros 



2006- 



Piran 



2004), are interesting as time-dependent candidate neutrino sources (jWaxman 



200 



i 



Meszaros k. Reeall993l ). Ice Cube conducts dedica te d searches wh i ch are tr iggered by satellite i n- 
formation for these objects (jAbbasi et al.l (j2010bl ). lAbbasi et al.l ()2009d ). (jAbbasi et al.l l2011al )). 
The untriggered all-sky search presented in this paper is also sensitive to this source class if two or 
more neutrinos can be detected from the same GRB. While the dedicated searches are in general 
much more sensitive (because they use the known time and direction of GRBs observed in gamma- 
or X-rays), the untriggered search has the potential to detect a burst which was not observed in 
photons (due to e.g. absorption or lack of monitoring). 

Another possibility for powerful emission is given by SGR's, X-ray pulsars that show vari- 
ability at different timescales and a persistent X-ray emission with luminosity L ~ 10'^^ erg/s 
with short bursts of X- and 7-rays with L ~ 10^^ erg/s lasting ~ 0.1 — 1 seconds (for review see 
Mereghettil (j2008l )). These X-ray pulsars, together with Anomalous X-ray Pulsars, are considered 
to be the best candidates for magnetars, isolated neutron stars powered by huge magnetic fields 
{B ~ 10^^ G). At times these sources emit giant flares with initial spikes of hard non-thermal 
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radiation up to luminosities of ~ 10^^ erg/s lasting some seconds. The enormous photon fluxes 
from t hese flares saturate most detectors. The flares may a lso accelerate baryons and produce neu- 
trinos (|Halzen et al-lbood : Eoka et al.|[2005l : IlIu et al.ll201ol ). Limits for photons in the 10 TeV-100 
PeV energy range using AMAN DA-II data were publishe d from the powerful giant flare observed 
in Dec. 2004 from SGR 1806-20 (|Achterberg et ahlboOGbl ). In the catalogue (Table H]) used in one 
of the triggered flare searches (Section [8]), we consider six days from a period of intense flaring 
from SGR 0501-1-4516 is included, which was discovered by SWI FT beginning Aug. 22, 2008 and 
also observed by RXTE/ASM, Konus-Wind and the Fermi GBM (jKumar et allboiol ). Even if the 
period of activity of SGR 0501-1-4516 is no t as intense as a giant flare and the observed power law 
component was quite soft (jRea et al.ll2009l ). this was a long-lasting period of flares that lead to the 
object's discovery as an SGR. Periods of activity from known SGRs IE 1514-57 and SGR 1806-20 
were omitted, since the large background of atmospheric muons severely limits IceCube sensitivity 
to southern hemisphere objects. 



3. The IceCube Detector and the Data Sample for this Analysis 



The IceCube Neutrino Observatory is composed of a deep array of 86 strings holding 5,160 digi- 
tal optical modules (DOMs), w hich are deployed between 1.45 and 2. 45 km below the glacial surface 
at the geographic South Pole (iHalzen Sz Hooper! 120021 : iBeckerl l2008l ) . IceCube strings are horizon- 
tally separated by about 125 m with DOMs positioned vertically 17 i n apart along each st ring. Each 



DOM consists of a Hamamat su photomultiplier w ith 25 cm diameter (jAbbasi et al.ll2010al ). electron 



tip. 

1 



ics for waveform digitization (jAbbasi et al.ll2009dl ). and a spherical, pressure-resist ant glass housing. 



IceCu be construction started with a first string installed in the 2005-6 season (jAchterberg et al 
2006al ) and has recently been completed in the austral summer of 2010-11. The configurations 



of IceCube that have been used for the present analyses (22-string and 40-string) are shown in 
Figure [H The IceCube Neutrino Observatory includes a dense subarr ay, DeepCore, designed to en- 
hance the physics performance of IceCube below 1 TeV (ICowenll2009l ) and a surface array, IceTop, 
for e xtensive air shower measurements on the composition and spectrum of cosmic rays (jStanev 



2010|). 



IceCube uses a simple multiplicity trigger, which requires that at least eight DOMs are triggered 
within 5 /US. For a DOM to trigger, it is both required that the DOM PMT voltage crosses the 
discriminator threshold (0.25 of a typical photoelectron), and this "hit" to be in coincidence with 
at least one other hit in the nearest or next-to-nearest neighboring DOMs on a string within ibl//s. 
This greatly reduces hits due to uncorrelated PMT noise and radioactivity in the glass. Once the 
simple multiplicity condition is satisfied, information from all triggered DOMs within a ±10 /is 
window is read out and merged to create an event. This means that 20 /is is the effective limit 
on how close two events can be in time for the 40-string or 22-string data. Improvements in event 
definition remove this constraint for data taken with the completed detector. Standard IceCube 
runs are eight hours long, runs in the 40 and 22-string configurations have a roughly two minute 
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Fig. 1. — : The growing IceCube detector seen from the top. Filled circles inside empty circles 
indicate deployed strings for each configuration, where all strings used in the 22-string configuration 
were also used in the 40-string configuration. 



gap between the end of one run and the beginning of the next, improvements have eliminated this 
gap in the operation of the full detector. The rate of each run is monitored a nd checked for an y 
deviation from an average that accounts for a seasonally adjusted average rate (iTilav et al.ll2009l ). 



Identification of neutrino-i nduced muon events in IceCube was demonstrated using atmospheric 
neutrinos as a calibration tool (jAchterberg et al.ll2007 ) . The measurement of atmospheric neutrinos 
by the 40-string configuration has been presented in (jAbbasi et al.ll2011bl ). This sample of upgoing 
events dominated by atmospheric neutrinos is used to look for astrophysical signals from point 
sources using directional and energy information. In the 40-string sam ple, the IceCube field of view 
has been extended compared to the upgoing only sample of 22-strings (lAbbasi et al.ll2009bl ) to also 
include downgoing events from the southern hemisphere. This technique was used for the firs t 
time to include events to —50° declination with the 22-string configuration (jAbbasi et al.l l2009al ) . 
The background sample in the downgoing region consists of very high energy muons. Atmospheric 
muons are roughly five orders of magnitude more numerous than neutrino-induced muons at the 
depth of IceCube. However, their number can be reduced by selecting high energy events so that 
the astrophysical signal can potentially emerge, if the signal spectrum is harder than that of the 
atmospheric background. This results in a different sensitivity for the northern hemisphere, where 
TeV-PeV neutrino astronomy is possible due to absorption of atmospheric muons in the Earth, 
with respect to the southern hemisphere, where only PeV-EeV neutrino searches are performed in 
the current analysis. 
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The selection of the data from the IceCube 40-str ing configurat i on use d in this search is 
discussed in the paper on the time-integrated searches, lAbbasi et al.l (|201ld ). We refer to this 
paper for all details on the selection of the muon events in the filtered stream (LI) that is sent 
over satellite from the South Pole and on the final cuts to obtain the data sample used in this 
analysis. The cuts were optimized using an E~'^ spectrum signal, as expected of neutrinos directly 
accelerated in astrophysical sou rces via stochastic shock accele ration. While th i s ide a was f irst 
presented by (jFermilll949l . Il954 ). and further developed by e.g. iKrymskiil ()l977l ): iBell (Il978al lbl). 
The resulting particle spectra follow a power-law close to E~'^. Detailed calculations show, however, 
that depending on t h e shock conditions , the spectra can also be somewhat flatter or steeper, see e.g. 
Stecker et al.l (j2007l ): iMeli et al.l (j2008l ). Here, we use an E~'^ spectrum as a first order estimate. 



It also includes a detailed discussion on the analysis method. The final sample consists of 36,900 
atmospheric neutrino and muon events from the whole sky (—85° to -|-85° in declination) detected 
by IceCube in the 40-string configuration in 375.5 days of good data taking, which corresponds 
to 92% uptime during the nominal operation period between April 5, 2008 and May 20, 2009 or 
Modified Juhan Date (MJD) 54561-54971. In this sample, 14,121 events are up-going, while 22,779 
events are down-going. Deadtime for this analysis is mainly due to test and calibration runs during 
and after the construction season. For time-dependent analyses in general, parts of the detector 
may be excluded for short periods from th e acquisition, but th e remaining part can be useful in 
case an astrophysical event occurs (see e.g. lAbbasi et al.l (j2009d )). 



In the triggered search for flares (Section [8]) we also consider events observed during data 
taking with 22 strings of IceCube, with a livetime of 275.7 days, or 89% of the ope ration period 



from M ay 31, 2007 to April 5, 2008 (MJD 54251-54561). This sample is described in lAbbasi et al 
(j2009bl ). and consists of 5114 candidate events from declinations -5° to -|-85°. 



(|201ld ) for the 40-string data and have also been discussed in lAbbasi et al. 



The systematic uncertainties have been evaluated and pre sented in Sec . 6 of Abbasi et al 

(|2009b|) for the 22- 



string data. The main uncertainties on the limits to the fluence of an E~'^ signal of muon neutrinos 
come from photon propagation, absolute DOM efficiency, and uncertainties in the Earth density 
profile and muon energy loss, accounting for a total of 16%. In this paper we focus on issues related 
to transient searches, such as the stability in time of the data sample and effects of the detector 
asymmetry for flares lasting less than one day. 

Before cuts are applied to the data, the samples are dominated by downgoing atmospheric 
muons. This is the case in the upgoing signal region as well, since some atmospheric muons are 
misreconstructed as upgoing and must be rejected in the process of applying analysis cuts. The 
atmospheric muon rate exhibits a seas onal variation of roughly ±10% due to changes in density 
of the atmosphere at the South Pole (jTilav et al.l |2009| ) . When the atmosphere is warmer and 
less dense during the austral summer, the fraction of pions and kaons in air showers which decay 
before interacting is increased compared to the winter. The muon rate also varies several percent on 
timescales of several days as a result of weather phenomena in Antarctica. For upgoing atmospheric 
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neutrinos the seasonal variations are smaller, approximately 5%, since neutrinos are created over a 
wide range of Earth's latitudes compared to the atmospheric muons created near the South Pole. To 
ensure stable detector conditions, the event rates of runs were required to be within 5cr from a rolling 
average over ±2 days. This loose constraint allows for short-term weather variability. All events 
have initial reconstructions performed using track and cascade based hypotheses, some of which 
are selected for transmission over satellite to the northern hemisphere for additional processing and 
analysis. The muon filter focuses on the selection of upward-going track-like events. The temporal 
variation of the event rate for the 22-string and 40-string runs is shown in Figure [21 where the 
seasonal modulation is clearly visible. The rate of the 40 string final event sample is shown in 
Figure [3] for upgoing and downgoing events and for the total sample. 
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Fig. 2. — : The rate per run of the filtered stream of muon events with zenith angle > 80° selected 
at the South Pole for the 22 (dark green points before MJD 54560) and 40 (gray points after MJD 
54560) string detectors as a function of MJD. The small modulations around the main seasonal 
oscillation are due to short-term weather variability (plotted in light blue). 

Due to the requirements for triggering and filtering, the cuts applied. Earth absorption proper- 
ties and detector geometry, the final sample of events is not uniform in the detector local coordinates 
zenith (0), and azimuth [(j)). For time-integrated point source searches, the azimuth dependence is 
usually neglected because it is smoothed in right ascension by the rotation of the Earth over long 
integration times. However, in a time-dependent analysis the azimuth dependence becomes impor- 
tant for time scales shorter than 1 day. The local coordinate (zenith and azimuth) distribution 
of 40 string data is shown in Figure [5] (left). In the northern sky there is the effect that events 
traveling along the longer axis (see Figure [T|) of the detector have a longer lever arm, and are more 
likely to trigger the detector and be well-reconstructed. Well reconstructed events typically have 
angular reconstruction errors of < 1° and energy reconstruction resolution of 0.3 in log(Energy). 
In the southern sky, there is an online cut on the integrated charge seen in all DOMs for a given 
event. This gives a preference to events which pass near a line of strings, yielding a six-fold peak 
in rates corresponding to the main axes of the detector symmetry. 
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Fig. 3. — : A graph of the rate of the final sample of 40 string events, in bins of 10 days. Errors 
are statistical. Also plotted are the individual rates of upgoing and downgoing events. The total 
fluctuation in the final data rate is ±5% for downgoing events and it ~ 4% for upgoing events. 
The analysis in Section [6] uses a sinusoidal fit to the atmospheric muon event rate to estimate rate 
fluctuations in the downgoing region. All other searches neglect these event rate fluctuations, which 
are negligible compared to statistical fluctuations in the expected signal. 



4. Unbinned Time-Dependent Likelihood Method 



The unbinned likelihood searches performed here are based o n the method descri bed in lBraun et al 



(j2008l ) and extended to searches for time-dependent behavior in lBraun et al.l (j2010l ). In this likeli- 
hood ratio method, a combination of signal and background populations is used to model the data. 
For a data set with N total events, the probability density of the i^^ event is given by: 



N 



Si + (l 



ns_ 

N 



)Bi. 



(1) 



where n<j is the unknown number of signal events with signal fraction of rig/N, Bi is the 
background probability density function (PDF) and Si is the signal PDF. The likelihood C of the 
data given the value of iig is defined as the product of the individual event probabilities: 



N 



^(-s) = n [^^^ + (1 - i)^^. 



(2) 



i=l 



This likelihood is maximized with respect to and any other fit parameters which are a part 



of the signal hypothesis. The maximization provides the best-fit values of these parameters. The 
background PDF, Bi, is given by: 



and is computed using the distribution of data itself. 

The spatial term B^^^'^^{6i,(j)i) is the event density per unit solid angle as a function of the 
local coordinates, shown in Figure H] (left). The energy probability, B™'^'^^^ [Ei^Oi), is determined 
from the reconstructed energy distribution of data as a function of the cosine of the zenith angle 



cos 6i ( see Figure U] on the right). This energy reconstruction, described in detail in lAbbasi et al 



(|201ld ). compares the measured to the expected density of photons along the muon track due to 
stochastic energy losses of pair production, bremsstrahlung and photonuclear interactions which 
dominate over ionization losses for muons above 1 TeV. The reconstructed energy value represents 
an estimate of the muon energy in GeV. The energy cut for the southern sky sample increases for 
smaller zenith angles, creating a strong zenith dependence of the energy in the southern sky as can 
be seen in Figured] (right). The goal of the energy cut is to sample a constant number of events 
per unit solid angle in the southern sky. Note that for the northern sky the energy dependence on 
zenith is small. The time probability Bl^'^^{ti,6i) of the background can be modeled including the 
expected seasonal modulations (done in Section [6]), which are less than ±10% and depend on the 
zenith angle, or taken to be flat since these modulations are negligible compared to possible signal 
fluctuations (all other searches). 

The signal PDF Si is given by: 



where S^^^^^ depends on the angular uncertainty of the event (Tj and the angular difference between 
the event coordinate Xi and the source coordinate Xg- S°^°^^^ , the energy dependent PDF which is 
a function of the reconstructed event energy Ei, and of the spectral index 7<j (a power-law spectrum 
with no cutoff such that dN/dE oc E~^'') is calculated from an energy distribution of simulated 
signal in a zenith band that contains the event. S\^^° is the time-dependent signal PDF. It depends 
on the particular signal hypothesis, which will be described in detail in each section. For each 
search, signal is injected with the same functional form in time as is being tested. Due to the 
low number of signal events involved, cross tests using one emission profile for the injection and 
a different one for the search yield largely similar results. The searches were designed to improve 
discovery potential with a minimum of source assumptions. 

The test statistic {TS) is calculated from the likelihood ratio of the background-only (null) 
hypothesis over the signal-plus-background hypothesis: 

TS = -2log\ = 1. (5) 
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Fig. 4. — : Left: The normalized event distribution of the 36,900 events in local coordinates for the 
40 strings data (the space term in Equation [3|) . There are two predominant effects: for upgoing 
events (northern sky, bottom half), events traveling down the longer end of the detector are more 
likely to trigger and pass cuts; for downgoing events (southern sky, top half), there are six peaks in 
the event rate. This is due to the initial filter conditions at the South Pole that select tracks more 
efficiently when they pass close to aligned strings. Right: The background energy distribution (the 
energy term in Equation [3]) , which is the normalized distribution of event energy values split into 
bins of the cosine of zenith. 

As expressed in Equation [Sj the test statistic values for scrambled samples will distribute as a 
chi-square function with number of degrees of freedom equal to the number of fit parameters. The 
best fit parameters hg , 7s and any fit time parameters Ts (which will be explained for each search) 
are obtained by maximizing TS. Several searches differ in the methods used to find the maximum 
value of TS, and may include weighting terms for values of a best fit parameter, which will be 
included in their description below. 

Larger values of TS are less compatible with the null hypothesis, and indicate its rejection 
at a significance level corresponding to the fraction of the scrambled trials above the TS value 
found in the data. Data scrambling is done by assigning a random time to each event from a 
period of active data taking and performing the proper coordinate transformation to get a new 
right ascension and declination. The fraction of trials above the TS value obtained from data is 
referred to as the p-value. This leads to the definition of the discovery potential at a particular 
significance level: the average number of signal events required to achieve a p-value less than a 
specified threshold in 50% of trials. IceCube uses a one-sided 5a significance level as the threshold 
for discovery, corresponding to a p-value of 2.87x10"'^. Similarly, the sensitivity is defined as 
the average signal required to obtain a p-value less than that of the median of the test statistic 
distribution of scrambled (background-only) samples in 90% of trials. 

Aside from the p-values from searches, in the absence of a signal upper limits on the fiuence 
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are provided, defined as the integral in energy and time of the flux upper limit: 



/ = / dt dEx E— = At dEx E^ = At^o [log(^max) - logl^^min)] , (6) 

where $o is the upper limit on the normalization on an E~'^ spectrum and Smin and Emax spec- 
ify the declination dependent energy range containing 90% of the signal spectrum at final sample 
selection criteria obtained from simulations. At is the duration of the emission. There is a cor- 
respondence between the fluence and the average number of events detected, shown as a function 
of the declination in Figure El The limits a re calculated according to the classical (frequentist) 



construction of upper limits (jNevmanI 119371 ) and the systematic error of 16% is neglected in all 



upper limits since the limits are dominated by statistical fluctuations for flares. 

Limits in this paper have been produced assuming a flux of only muon neutrin os. The scenario 



of standard neutrino oscillations over astronomical distances (jAthar et al.l 120001 ) assumes equal 
fluxes of all flavors of neutrinos at the Earth from a source producing neutrinos via pion decay 
with a ratio of z^e : : ^^r = 1 : 2 : 0. When considering equal fluxes of muon and tau neutrinos 
at the Earth, the resulting upper limits on the sum of both fluxes are about a factor of 1.7 times 
higher than if only muon neutrinos are considered rather than the expected factor of two due to 
oscillation if no tau neutrinos would be detectable. For an E~'^ spectrum of the signal neutrino flux 
the contribution due to the detectable tau neutrino flux for sources at the horizon is 10% and up to 
15% for sources in the Northern hemisphere. This is due to the tau decay channel into muons with 
a branching ratio of 17.7% and in part to the tau leptons with energy greater than some PeVs that 
may travel far enough to be reconstructed as tracks in IceCube before decaying. In the upgoing 
region we have considered tau regeneration in the Earth. 



5. All-Sky Time Scan 

The all-sky time-dependen t search presented h ere complements the all-sky search applied to 



the IceCube 40-string data in (jAbbasi et al.l l201ld ). While that search has the best sensitivity 



to steady sources, a source which has emitted neutrinos for only a limited period of time might 
not be detected. The time-dependent analysis here scans for a significant excess with respect to 
background over all time scales (from sub-seconds to the full year) at each direction of the sky. 
For flares shorter than '^lOO days, the discovery potential of the time-dependent search typically 
becomes better than the time-integrated one, and in principle a short burst can be discovered 
with only two events if they occur close enough together in time (~ 0.1 seconds for E~'^ spectrum 
events). The advantage of such untriggered searches is their ability to cover all emission scenarios, 
including neutrino emission without any observed counterpart in the electromagnetic spectrum. 
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Fig. 5. — : The total fluence from an E~'^ spectrum muon neutrino signal in a declination band 
divided by the number of events in the band in the 40-string configuration, plotted against decli- 
nation. 



5.1. Method and Expected Performance 



The method in lBraun et al.l ( 20ld ) is adapted for this search to a real detector with non- uniform 
acceptance and deadtime. The non-uniform acceptance can be seen in IH deadtime compensation 
is shown below in O The time-dependent probability density function from Equation 3] for this 
search is a Gaussian function: 



where ti is the arrival time of the event, and fit parameters To and aj- are the mean and sigma of 
the Gaussian describing flaring behavior in time. The maximization of the test statistic returns the 
best-fit values of the Gaussian mean (the time at which the flare peaks) and sigma (corresponding 
to the duration of the flare). Both the background and expected number of events are small, 
distinguishing a box-type function from a Gaussian would require many more events than required 
for a 5f7 discovery, and we find that using either of these flare hypotheses performs similarly (see 
Section E]). It was found that the fitting method used in this section worked better with a continuous 
function, so a Gaussian functional form was chosen. 



Because there are many more independent small time windows than large ones, the test statistic 
formula of Equation [5] is modified to include a weighting term to correct fo r this effective trial factor 
and avoid undue preference for short fiares using a Bayesian approach (jBraun et al.1120101 ). The 
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test statistic formula that is maximized is then: 



TS = -2 log X ^("- = °l 1 , (8) 

where the first factor in the square brackets is the weighting term and the second is the likelihood 
ratio. T is the total livetime of data taking, ns,7, (Tt,To are the best-fit values for the number 
of signal events, spectral index, width and mean of the Gaussian flare, respectively. In order to 
prevent the weighting term from becoming less than 1, a constraint is placed on the flare width 
aT- This is done to prevent flares with zero amplitude (n<j=0) from having a positive test statistic, 
which would happen if the flare width ax were allowed to be greater than T/\/27r. 



As described in lBraun et al.l ([20101), the numerical maximizer needs an initial candidate flare 
(a "first guess"). In that paper this first guess was obtained by selecting events within 5° of the 
source location (in this analysis we use the criteria Si/Bi > 1, where Si and Bi are defined in 
Section U omitting the time term). The data is broken up into sets of m temporally consecutive 
events, where 2 < ?n < 5, for initial testing. So, for a stream consisting of time-ordered events 
numbered 1,2,3,4,5,6,7..., the initial test uses events (1,2), (2,3), (3,4) etc., (1,2,3), (2,3,4), (3,4,5) 
etc., (1,2,3,4), (2,3,4,5), (3,4,5,6) etc., and (1,2,3,4,5), (2,3,4,5,6) etc. Each set is tested using the 
described likelihood formula for compatibility with a flare with an spectrum. The candidate 
with the best test statistic (from Equation [8]) is used as the initial first guess in the maximization. 
In the current analysis, t he maximurn numb er of consecutive events for the initial test has been 



increased with respect to lBraun et al.l (j2010l ) to m = 10, improving the sensitivity to longer flares. 
This brings the performance of the analysis close to that of the corresponding time-integrated 
analysis at large time scales. Given that more than 5 events are required for discovery for > 2 
days (see Fig [7]) , if the maximum is not increased the method will occasionally only find a subset of 
the injected events, hence increasing the total signal required to cross the threshold for discovery. 

Figure [7] (left) shows the mean number of injected events from a Gaussian time function needed 
for a 5(7 discovery for 50% cases (black solid line) as a function of the duration of the flare cjt for 
a fixed source location at declination of +16°, though sources at other declinations yield similar 
results. This is compared to the number of events needed in a time-integrated search (black dashed 
line): the number of events needed to discover a flare of Is duration is about a factor of 4 lower 
than for a time-integrated search. At long timescales the flare search performs only 10% worse 
than the time-integrated search, even with 2 additional free parameters in the fit. In the same plot 
the median upper limits at 90% c.l. are shown for the time-dependent search and for the time- 
integrated one. On the right the corresponding fluence is given, where a correction is introduced 
for the median dead time during a given flare as a function of the flare width (see Figure [6|) . 

The fact that the 50% 5a discovery potential curve descends below the 90% median upper 
limit curve is due to the effect of Poisson statistics. The untriggered search must observe at least 
two events in order to identify a flare. For a simulated flaring source which injects a mean number 
of events fi, fi must equal at least 1.68 for 50% of simulated trials corresponding to 2 or more signal 
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Fig. 6. — : The fractional duration of randomly-simulated flares which occur during the uptime 
for the 40-string configuration for a range of different flare durations. The black line marks the 
median fraction of fluence occurring during the detector livetime for a given flare duration, which 
is used as a correction factor for the fluence of observed flares. For instance, for flares shorter than 
one minute, there is approximately an 8% chance of the flare occurring completely during detector 
downtime. Flares longer than one day will always have some emission during uptime; on average 
92% of the total emission will coincide with usable run time. 
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Fig. 7. — : The 50% 5a discovery potential and 90% median sensitivity in terms of the mean number 
of events (left) and fluence (right) for a fixed source at +16° declination. The number of events 
for the median sensitivity and discovery potential for the time- integrated search are also shown. 
Flares with a or of less than 100 days, or a FWHM of less than roughly half the total livetime, 
have a better discovery potential than the steady search. 



events. Therefore, at the shortest timescales, the mean signal needed for a discovery in 50% of trials 
asymptotically approaches 1.68 events. We find the sensitivity at 90% CL saturates at 2.9 events, 
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which is aheady near the time-independent sensitivity of 3.15 events and the statistical hmit. This 
is the reason why the discovery potential curve is lower than the sensitivity in Figure [71 

The method is applied as an all-sky scan over a grid (0.5° x 0.5°) in right ascension and 
declination. The final result of the analysis is the set of best fit parameters from the location with 
the highest test statistic value. A final value for this analysis is obtained by performing the same 
scan on scrambled data sets, and counting the fraction of scrambled sets which have a maximum 
test statistic greater than or equal to the maximum found in the data. 



5.2. Results 

Using the 40-string data, the location which deviates most from the background expectation is 
found at (RA,Dec)=(254.75° , -|-36.25°). Two events are found, with a best-fit spectrum 7 of 2.15, 
mean of the fiare of MJD 54874.703125 and width (fx of 15 seconds. The two events are 2.0° 
apart in space and 22 seconds apart in time. The — logio(^'- value) corresponding to this observation 
is 4.67. A clustering of higher significance is seen in 56% of scrambled skymaps (276 out of 500), a 
result consistent with the null hypothesis of background-only data. 

Figures [8l to 1 101 show maps of the pre-trial p- values and best-fit parameters To and cfr- Figures[9] 
and [10] require that the best-fit number of signal events be greater than zero, white area corresponds 
to being consistent with no flare being detected. A large effective trial factor of 2.6 x 10^ is generated 
by scanning the whole sky. Therefore it is desirable to look only at a few sources in order to decrease 
the trials, which is done in Section [6l 




Fig. 8. — : The equatorial coordinate map shows the p-value of the most significant fiare in time 
and space at each location of the grid where the likelihood is calculated. The p-value is indicated 
on the z-scale on the right. The black curve is the Galactic plane. 



- 20 - 



6. Time Scan for Candidate Sources 



By targeting specific, a priori promising directions in the sky, an analysis can reduce the 
effective trial factor of the all-sky scan. One way this can be done is by performing the analysis 
in Section [5] at the specified locations only. Here, thi s sear ch wa s instead implernented using a 
time-clustering algorithm developed in ISatalecka et al.l (120071 ) and iBazo-Alba et al.l (j2009l ). which 
achieves similar performance. The algorithm finds the most significant flare in a period by testing 
the most promising time windows, which are defined by the times of the neutrino events. 

For the source list, variable bright astrophysical obj ects from the enti re sky are selected. 
Sources are taken from the Fermi-LAT Bright Source List (jAbdo et al.ll2009bl ). whose data taking 
period overlaps with the IceCube 40-string sample. The sources include 30 blazars (24 FSRQs 
and 6 BL Lacs), one high- mass X-ray binary, one radio galaxy and seven unidentified objects. In 
this analysis the following selection criteria are defined for choosing the most promising variable 
astrophysical sources: 



Classified as variable by Fermi-LAT, 

Flux [100 MeV - 1 GeV] > 1.1 x lO"'^ photons cm'^s"!. 



The definition of variability provided by the Fermi-LAT Bright Source list is that the obser- 
vation has a probability of less than 1% of being a steady source (i.e. variability flag=T). The 
second requirement sets a minimum photon flux, motivated by the correlation between neutrinos 
and photons emitted from the source predicted by hadronic models. The average photon flux from 
100 MeV to 1 GeV of all Fermi variable sources is 2.3 xlO~^ photons cm~^s~^. The flux threshold 
chosen keeps 60% of these sources. The list of selected candidates contains 40 objects (see Tabled]), 
18 in the southern sky and 22 in the northern sky. 



6.1. Method and Expected Performance 

For a given source location, signal- like events are defined as having a time-integrated Si/Bi > 1, 
where Si and Bi are defined in Section [U omitting the time term. Each pair of these event times 
assigns a starting and ending time (ij and tj, respectively), to the fiare search windows ATy = tj—ti. 
The longest flare duration is constrained in the algorithm to be 30 days. Apart from this constraint 
the algorithm loops over the events, testing all windows deflned by each set of signal-like events 
specified above, maximizing for the signal fraction and spectrum for each set. The test statistic is 
calculated for the most promising fiare time windows and the best time window (i.e. the highest 
test statistic value) is chosen as a flare candidate. 

The signal time probability 5**"^° is deflned by: 
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Sf^^{U,t,) = -^ . (9) 

The time probability 5'*™'=^ is constant, no particular time structure is assumed during the flare. 
The statistics of signals expected will be small enough that no particular functional form should 
be discernable. 

The mean number of events needed for this analysis to achieve a Scr discovery with 50% of 
trials was calculated for different widths of simulated flares (see Figure [TT]) . The flare duration was 
investigated in the range from 30 days to 20^s, the minimum time between events. The discovery 
potential is very similar to the method of Section IH which is also plotted in Figure [TTl For flares 
with duration on the order of minutes, one third of the events needed in a time integrated search 
is necessary for a detection with the untriggered flare method. 



6.2. Results 

The time-scan looking for neutrino flares was applied to the 40 selected source candidates 
using IceCube 40-string data. No signiflcant excess above the atmospheric background is found. 
The results and upper limits for each source are presented in Table [T] and summarized in Figure 
[HI The highest fluctuation observed corresponds to OFGL J0643.2+0858 (dec=8.9°, ra=100.8") 
with a p-value of 7% (1.5 a) pre-trial. The corresponding best time cluster was 14.3 days, lasting 
from MJD 54846.5 to 54860.8. Correcting for the trial factor from looking at 40 sources, the final 
post-trial p-value is 95%. The post-trial p-value is obtained from the distribution of the maximum 
test statistic for many equivalent samples obtained by scrambling the time of events. 
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Fig. 9. — : The equatorial coordinate map shows the best fit of the mean time of the flare Tq 
(MJD-54,000) for the most significant flare found at each location of the grid where the likelihood 
is calculated. The black curve is the Galactic plane. 




Fig. 10. — : The map in equatorial coordinates of the best fit width cj', in days, of most significant 
flare at a given location found at each location of the grid where the likelihood is calculated in the 
search. The black curve is the Galactic plane. 



- 23 - 



03 
> 

o 



10 



03 

E 



CO 
03 




10' 



Time-integrated searcti 
Time-Clustering Algorithm 
Gaussian Time-Dependent Search 



10' 



10-^ 



^ Width of flare (days) 



Fig. 11. — : Mean number of events for a 5cj discovery in 50% of trials as a function of the flare 
duration, calculated as an example for a point source at dec=16° ra=343° using the time-clustering 
method (this section) and the method from section [H 
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Fig. 12. — : Fluence upper limits (red triangles) of the most significant cluster for each of the 
40 selected sources calculated in the time windows given in Table [1] versus declination. The blue 
squares represent the median sensitivity on the fluence calculated for the same time windows. 



Table 1:: Results for pre-defined variable astrophysical source candidates using the likelihood time clustering algorithm. 
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675 


J0643.2+0858 




8.983 


100.8 


3.7 


6.7 


0.07 


14.3 


54846.5 


2.57 


J2147.1+0931 


PKS 2144+092 


9.519 


326.8 


3.7 


6.6 


0.49 


27.8 


54737.9 


1.1 


J1751.5+0935 


OT 081 


9.591 


267.9 


3.7 


6.6 


0.66 


6.22 


54917.2 


0.84 


J2327.3+0947 


PKS 2325+093 


9.794 


351.8 


3.7 


6.6 


0.82 


2.84 


54603.2 


0.339 


J1504.4+1030 


PKS 1502+106 


10.51 


226.1 


3.7 


6.6 


0.17 


6.47 


54777.4 


1.32 


J1553.4+1255 


PKS 1551+130 


12.92 


238.4 


3.7 


6.5 


0.48 


0.488 


54566.3 


0.885 



Table 1 — continued from previous page 



Source 


Other Name 


dec [°] 


ra n 


log Erain 


log ^max 


value 


At 


^start 


Fluence Limit 


(OFGL) 








(GeV) 


(GeV) 




(days) 


(MJD) 


(GeV/cm2) 


J0531.0+1331 


PKS 0528+134 


13.53 


82.76 


3.6 


6.4 


0.38 


0.581 


54597 


1.02 


J2254.0+1609 


3C 454.3 


16.15 


343.5 


3.6 


6.3 


0.66 


7.94 


54594.4 


0.71 


J0238.6+1636 


AO 0235+164 


16.61 


39.66 


3.6 


6.3 


0.12 


0.216 


54776.6 


1.25 


J1522. 2+3143 


TXS 1520+319 


31.73 


230.6 


3.4 


6.0 


0.52 


3.95 


54869 


2.11 


J1310.6+3220 


B2 1308+32 


32.34 


197.7 


3.4 


5.9 


0.76 


27.7 


54671.8 


0.716 


J1635.2+3809 


4C +38.41 


38.16 


248.8 


3.3 


5.8 


0.092 


0.268 


54795 


1.66 


J1641.4+3939 




39.67 


250.4 


3.3 


5.8 


0.29 


0.268 


54795 


1.5 


J0320.0+4131 


NGC 1275 


41.52 


50 


3.3 


5.8 


0.62 


28.3 


54576.7 


1.26 


J0222. 6+4302 


3C 66A 


43.04 


35.65 


3.3 


5.7 


0.51 


19.6 


54641.1 


1.4 


J0654.3+4513 


B3 0650+453 


45.22 


103.6 


3.3 


5.7 


0.93 


5.85 


54903.5 


1.69 


J0240.3+6113 




61.23 


40.09 


3.1 


5.5 


0.91 


6.19 


54699.5 


1.71 


J1849.4+6706 


S4 1849+67 


67.1 


282.4 


3.1 


5.3 


0.7 


10.6 


54708.9 


1.31 


J0722.0+7120 


S5 0716+71 


71.35 


110.5 


3.0 


5.3 


0.75 


7.11 


54864.7 


1.47 



Note.- The source name is the OFGL catalogue designation. The p-value was calculated from simulated background skymaps, At is the 
flare duration of the best cluster and Tgtart its starting time. The fluence upper limit was calculated by integrating d^/dE x E over the 
declination dependent energy range from i?min to -Emax to contain 90% of signal spectrum and At, assuming a neutrino energy spectrum of 

E-\ 
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7. Triggered Searches Based on Continuous Photon Observations 



When there is specific timing information about the activity of an astronomical object, that 
information can be used to perform a targeted search with reduced background. This section 
describes searches in which the photon observations are essentially continuous, and this complete 
set of flux measurements in time is used. For flares lasting on the order of one day, MWL information 
can produce a discov ery with about one third fewer signal events with respect to untriggered searches 
teraun et all (l20ld )). 



The source selection was motivated by Fermi alerts, which are only issued for sources seen at 
a flux level greater than 2 x 10^^ photons cm~^. The selected sources are listed in Table [2j 
These sources include 6 FSRQs, one BL Lac and one unidentified object. The lightcurves were 
produced for this work using the Fermi Public Release data, using the diffuse class event selection. 
For each source the Fermi Science Tools v9rl5p2 package is used to select photons from within 2° 
of each source and calculate the total exposure. Photon events with zenith angles greater than 105° 
were excluded to avoid contamination due to the Earth's albedo. Time bins of one day width were 
then used to calculate an average flux. There are two modifications to thi s procedure : the b lazar 
3C 454.3 was seen in a massive outburst before official science operation (jAbdo et al.ll2009al ). for 
this source the published lightcurve is used. Also, the source PKS 1502+106 was noted to have 
a large outburst immediately before official s cience operations beg an, extending several days after 
the public information begins (|Ciprinill2008l : lAbdo &: others! l2010l ). PKS 1502+106 is taken to be 
flaring since the time of the alert at a fixed flux level. This flaring activity is a possible orphan flare 
in hard X-rays, since the SWIFT-BAT did not observe any evident flare in the 15-50 keV band 
while SWIFT XRT and UVOT observed a flare in soft X-rays and optical. 



7.1. Method and Expected Performance 



A Maximum Likelihood Block (MLB) algorithm (IScargle 



1998 



Resconi et al.ll2009l ) is used to 



denoise the lightcurves by iterating over the data points to select periods from the lightcurves which 
are consistent with constant flux once statistical errors are taken into account. The MLB algorithm 
compares the likelihood that a set of points between Xi and xj is compatible at a confidence level 
with the change of state between the points Xc and Xc+i contained in the interval. The confidence 
level requires that for a given set of data points from Xi to xj that the set of points from Xi and Xc 
and Xc+i to Xf be: 



C{Xi,Xc)C{Xc+l,Xf) 

log — ^ — > log C 



(10) 



C{xi,Xf) 

where C{xa,Xb) represents the likelihood that a set of points from Xa to x^ represents a constant 
flux state from the source between the points, and C is the confldence level. The level of the flux 
state is determined using the error-weighted mean of the points tested. The method iterates over 



-27- 



the different possible changepoints Xc, taking the most likely changepoint for the entire dataset and 
iterating over each subsection of the data. We tested values of logC from 1 to 1000, for values 
below 5, the typical denoised lightcurve typically follow each data point, for values from 9 to 100 
very similar results for the denoised lightcurves were found. The final value of log C used in the 
analysis was 20. 

With the hypothesis that the intensity of the neutrino emission follows the intensity of the 
photon lightcurve, the signal time PDF is simply the normalized lightcurve itself. A slightly modi- 
fied hypothesis is that the neutrino emission follows the lightcurve, but only when the photon flux 
goes above a certain threshold i<th- In this case, the value of -Fth can be used as a free parameter in 
the analysis, finding the value of the threshold which maximizes the significance of the data. This 
method also avoids any penalty from making an incorrect a priori choice on a flaring threshold. 
F{ti) is defined as the value of the denoised lightcurve at ti and the fit parameter Fth is the flux 
threshold below which no neutrino emission is assumed (i.e. 5'|™'^=0 if F{ti) < Fth)- In the case 
of F{ti) > Fth, the probability of neutrino emission is assumed to be proportional to the flux level 
above that threshold: 
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Fig. 13. — : (Left) An example of the one-day binned Fermi lightcurve (blue points, with statistical 
errors) and denoised lightcurve (pink solid line) for the blazar PKS 1510-089. The dashed line 
is an example fit threshold. The lightcurve begins here on August 10, 2008 (MJD 54688), when 
Fermi science operations began, while the time axis shows the entire 40-string data taking period. 
(Right) The time PDF used in the neutrino signal hypothesis corresponding to the example photon 
threshold shown in the left graph (5 x 10~^ photons s~^ cm~^). 



- N} ' ^^^^ 
where the normalization factor Nj is the integral of the denoised lightcurve above the threshold. 
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This time-dependent PDF is then used as before in Equation [H This method is iUustrated in 
Figure [ni 
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Fig. 14. — : The plot of the 5cj 50% discovery potential for the source PKS 1510-089 (the corre- 
sponding hghtcurve is shown in Figure [T3|). as a function of the true flux threshold for neutrino 
emission (left) and as a function of the duration the lightcurve spends above the threshold (right). 
The discovery potential curves are plotted for the time- integrated case (black short dashed line), 
and from bottom to top for the case where the threshold is fixed to the true threshold (solid red 
line), the case where the threshold is a free parameter (black long dashed line, used in this analysis) 
and the case where there is an unknown lag (up to it 50 days) between GeV and neutrino emission 
(blue dashed line). 

The effect of adding this additional degree of freedom by fitting for F^^^ can be seen in Figure [TTl 
The effect is small compared to the penalty of fixing the threshold to an incorrect value. The effect 
of allowing an unknown lag up to ±50 days between the photon and neutrino emissions was also 
tested, and was found to give a marked increase in the number of events required for discovery. 
Hence, we used the method allowing only up to a ±0.5 day lag that accounts for the 1 day binning 
of lightcurves. 



7.2. Results 

The results from all sources are listed in Table[2l The most significant source is PKS 1502+106, 
which has a pre-trial p- value of 5%. The method finds one high-energy event during the August 
2008 flare. The prescription to provide the post-trial p- value was to consider the most significant 
among the flares obtained in this and Section [8] with the 40-string configuration. The post-trial 
p- value is 29%, which is compatible with background fluctuations. 
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8. Triggered Searches Based on Intermittent Photon Observations 

Ground based observatories such as H.E.S.S., MAGIC, and VERITAS cannot monitor sources 
continuously, because they can only operate when there is good nighttime visibility. Their obser- 
vations are nevertheless extremely important for neutrino searches, because they detect photons at 
TeV energies that are potentially better correlated to neutrinos of the energies to which IceCube is 
sensitive. While these observatories can issue alerts for source activity, they often cannot constrain 
the beginning or end of the activity to within a few days. For alerts such as these, the present 
analysis uses a simple time cut, taking a window for events one day before and after the identified 
flare. The catalogue corresponding to the 40-string data includes a recent s uspect "orphan flare" 



at the level of 10 C rab from Mrk 421 observed by VERITAS and MAGIC (jAcciari et al.ll2009al : 



Vittorini et al.ll2009l ). During the 40-string period, there was no overlap between sources with flares 



seen by Fermi and by lACTs and tested in this paper. 

8.1. Method and Expected Performance 

The nature of this analysis is a simple cut in time between tmin and tmax, which can be expressed 

as: 

^time _ -^(^max ~ U) ^ H{ti — tmin) q2'\ 

* ~ f — f ■ 

''max ■'mm 

where ti is the arrival time of the event, tmax and tmin are the predefined the upper and lower 
bounds of the time window defining the flare, and H is the Heavyside step function. This time- 
dependent signal PDF is then used in Equation [H In this analysis, the signal population size 
and spectrum index 7 are the only fit parameters. 



8.2. Results 

Five of the seven flares tested with the 22 string data (Table ED showed no excess of events in 
the vicinity of the corresponding sources in the selected time windows, while S5 0716+71 and lES 
1959+650 showed one event each. The post-trial p- value is 14%, the most significant fiare being 
the 10 day flare of S5 0716+71. This result is consistent with background fluctuations. 

Of the six sources tested with the 40-string data (Table U]), flve showed no excess of events in 
the vicinity of the sources during the selected time periods. The flnal post-trial p-value for the 40 
string analysis (considering the 15 flares in Tableland Table |4| is 29%. 
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9. Conclusions 

In this paper we discuss four time-dependent searches: two are "untriggered" and scan over 
direction, energy and time to find clusters of neutrino events; the others are "triggered" by multi- 
wavelength information. While the first two searches are generic and sensitive to flares not seen 
in photons, the others are more sensitive because of the reduced trial factor but concern specific 
catalogues of variable sources. Time-dependent searches can be more sensitive to short flares thanks 
to the reduction of the background of atmospheric muons and neutrinos over short time scales. The 
untriggered search using a predefined catalogue of 40 variable sources shows a post-trial p- value of 
95% and upper limits on the fluence are calculated (see Table [1]). Eight of these sources are also 
triggered on observed flares in the triggered search. The all-sky scan over all directions finds that 
the most significant cluster of events is separated in time by 22 s and in space by 2° and has a 
p-value of 56%. The most significant observation of fiaring from 14 sources in catalogues compiled 
using Fermi-LAT and lACT alerts during the 40-string configuration data taking is PKS 1502+106, 
with a p- value of 29% after trials. The most significant fiare triggered by MWL observations during 
the 22-string configuration is S5 0716+71 with a p-value of 14% after trials. All these results are 
compatible with a fluctuation of the background. 

The complete IceCube detector began taking data in April 2011, and is expected to have an 
effective area twice that of 40 strings at high energies, and up to a factor of ten at 100 GeV. This 
is especially a boon to time-dependent neutrino point source analyses, which are limited by the 
statistics of the signal events observed. Time-dependent analyses using data from the full detector 
will have roughly a factor of two to three improvement in the fluence upper limits and discovery 
potentials from point sources, a significant improvement in the capabilities of neutrino astronomy. 
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Source 


dec [°] 


ra [°] 


log Emin 


log E'max 


pre-trial 


Best-Fit Threshold 


Duration 


Upper limit 








(GeV) 


(GeV) 


p- value 


(10^^ cm"2 s-i) 


(days) 


(GeV/cm2) 


PKS 1510-089 


-9.1 


228.2 


4.8 


7.8 







282 


2.49 


3C 66A/B 


43.0 


35.7 


3.3 


5.7 


0.47 


0.675 


57 


0.778 


3C 454.3 


16.2 


343.5 


3.6 


6.3 


0.20 


9.47 


2.5 


0.422 


PKS 1454-354 


-35.6 


224.4 


5.6 


7.8 







282 


13.1 


3C 279 


-5.8 


194.0 


4.2 


7.6 


0.47 


2.34 


6 


1.69 


PKS 0454-234 


-23.4 


74.3 


5.5 


7.9 







282 


7.87 


PKS 1502-^106 


10.5 


226.1 


3.7 


6.6 


0.049 


3.13 


8 


0.370 


J123939-h044409 


4.7 


189.9 


3.7 


7.1 







282 


0.661 



Table 2:: Sources tested with the 40 string data and pre-trial p-values for the flare search with continuous lightcurves. In the 
event of an under-fluctuation no p-value is calculated and there is no indication of flaring behavior. The overlap between the 
Fermi public release data and the 40-string data taking period is 282 days. The duration column corresponds to the amount of 
time for the lightcurve which is above the best-fit threshold. 



Source 



dec [°] 



ra 



Alert Ref. 



Time 
Window 
(MJD) 



log E^min 
(GeV) 



log E^max 

(GeV) 



p- value 
(pre-trial) 



lES 1959+650 65.1 300.1 

lES 1959+650 65.1 300.1 

3C 454.3 16.2 343.5 

3C 454.3 16.2 343.5 

Cyg X-1 35.2 299.5 

S5 0716+71 71.3 110.5 

S5 0716+71 71.3 110.5 



( Bottacini. E. et al.ll2007l ') 



( WhipplellioOTh 



( Vercellone et al.ll2008h 



( Vercellone et al.ll2009h 



( Golenetskii. S. et al.ll2007l ) 



( Chen. A. W. et al.lliooi ') 
( Chen. A. W. et al.lliooj ') 



54428- 
54435.5- 
54305- 
54416- 
54319.5- 
54350- 
54392- 



54433 

54440.5 

54311 

54426 

54320.5 

54356 

54402 



3.1 
3.1 
3.8 
3.8 
3.4 
3.1 
3.1 



5.4 
5.4 
6.5 
6.5 
5.0 
5.3 
5.3 



0.08 



0.02 



Table 3:: Flare list during the 22 strings data-taking: source name, references for the alert, interval in modified Julian day, 
pre-trial p-value. The p-value is calculated only when hs is greater than zero. 



Source 


dec [°] 


ra n 


Alert Ref. 


Time 


log Emin 


log E^max 


p- value 


Fluence 










Window 


(GeV) 


(GeV) 


(pre-trial) 


Limit 










(MJD) 








(GeV/cm2) 


Markarian 421 


38.2 


166.1 


fPichel 2009^ 


54586-54592 


3.3 


5.8 




1.51 








fVittorini et al. 2009) 


54621-54631 










W Comae 


28.2 


185.4 


(Acciari et al. 2009bl 


54623-54627 


3.4 


6.0 




1.32 


S5 0716+714 


71.3 


110.5 


fMazin et al. 2009^ 


54572-54582 


3.0 


5.3 


0.34 


3.26 


SGR 0501+4516 


45.3 


75.3 


fPalmer 2008) 


54700-54706 


3.2 


5.7 




1.72 


lES 1218+304 


30.2 


185.3 


fAcciari et al. 2010) 


54859-54864 


3.4 


6.0 




1.43 


Markarian 501 


39.8 


253.5 


fPichel 20091 


54951-54953 


3.3 


5.8 




1.78 



Table 4:: Flare list seen with occasional coverage during the 40-string data-taking. References are for the alert which prompted 
the selection. The p-value is calculated only when hg is greater than zero. The flare windows for Markarian 421 were added 
together, only one p-value and upper limit for both periods is calculated. 



